Experimentally accessible reentrant phase transitions in double-well optical lattices 
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We study the quantum phases of bosons confined in a combined potential of a one-dimensional 
double-well optical lattice and a parabolic trap. We apply the time-evolving block decimation 
method to the corresponding two-legged Bose-Hubbard model. In the absence of a parabolic trap, 
the system of bosons in the double- well optical lattice exhibits a reentrant quantum phase transition 
between Mott insulator and superfluid phases at unit filling as the tilt of the double- wells is increased. 
We show that the reentrant phase transition occurs also in the presence of a parabolic trap and 
suggest that it can be detected in experiments by measuring the matter-wave interference pattern. 
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Recently, the system of ultracold bosonic atoms in 
a double-well optical lattice (DWOL) has been used 
to study quantum information processing and quantum 
many-body systems 0, H S 0] . A DW0L 

is constructed 

from two-color two-dimensional (2D) lattices in the hor- 
izontal xy plane and another lattice in the vertical z di- 
rection, resulting in a 3D lattice whose unit cell is a dou- 
ble well. By exploiting precise control of the properties 
of the double wells, recent experiments have realized an 
atom interferometer (2j, two-atom exchange oscillations, 
which were used to create a quantum SWAP gate [|[ , and 
correlated tunneling of two interacting atoms Q ■ 

In single-well optical lattices, the transition from a 
superfluid (SF) phase to a Mott insulator (MI) phase 
has been observed by adiabatically increasing the lat- 
tice depth [1, @]. In previous DWOL experiments all 
the double wells were independent of each other due to 
the large potential barrier between them, and the sys- 
tem was always in an insulating phase. However, the 
potential barrier could be lowered to induce a transition 
to a SF phase In particular, one can investigate the 
SF-MI transition in a two-legged ladder ID DWOL by 
tuning the lattice depth in the vertical direction. While 
in single- well lattices the SF-MI transition is determined 
by the ratio, U/t.oi the on-site interaction energy to the 
hopping energy @, H[ , in DWOL the SF-MI transition is 
governed by the competition between the on-site interac- 
tion and the tilt (energy offset) of the double wells. 

In this paper, we study the ground state properties 
of bosons in a two-legged ladder lattice with a parabolic 
confining potential. This case is more representative of 
current experiments than is an unconfined, homogeneous 
lattice. Although mean-field-type methods describe the 
SF-MI transition of bosons in optical lattices very well 
in 3D [7[, they fail to describe quantitatively the two- 
legged ladder system due to strong quantum fluctua- 
tions 8]. Hence, we use the time-evolving block deci- 
mation (TEBD) method [1G3, 

which is one of the best 
methods currently available to study strongly correlated 
ID quantum systems. We show here that a reentrant 



FIG. 1: (color online) Schematic double-well lattice potential 
in the xz plane for a = 0.2, = 0, showing the existence of 
isolated two-legged ladders. 



phase transition (MI-SF-MI) occurs as the tilt of the dou- 
ble wells is changed for a fixed parabolic confinement. 
Moreover, a distinctive signature of the reentrant phase 
transition appears in the momentum distribution, which 
is directly related to matter-wave interference patterns 
observed in experiments. 

We consider a Bose gas at zero temperature confined 
in a combined potential of a DWOL and a parabolic trap. 
A DWOL potential in Refs. [J 0, Q is given by 

V(r) = V in (sin 2 (fez) + sin 2 (ky) - 2) 

-V out (sin(fcc - 0) + cos(fcy)) 2 + V z sin 2 (fcz), (1) 

where Vi n is the depth of the "in-plane" lattice with pe- 
riod of d = n/k, whose light polarization is in the xy 
plane, V out the depth of the "out-of-plane" lattice with 
2d-periodicity [l[ , and V z is the depth of the lattice in the 
vertical z direction. When the ratio a = V^ ut /Vj n satis- 
fies < a < 0.5, the unit cell of the potential Eq. ((T|) 
is a double-well. Control of the ratio a and the relative 
phase 8 provides the flexibility to adjust the double- well 
parameters: the barrier height and the tilt. When = 0, 
the double wells are symmetric (no tilt). Increasing a 
from zero enlarges the barrier height Vhigh between dou- 
ble wells in the xy plane and reduces the barrier height 
Mow between two wells within a double well. Keeping 
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Vhigh — Mow ^ En and tuning T4, one can create an en- 
semble of ID DWOLs with a two-legged ladder structure, 
where all the two-legged ladders are decoupled from each 
other as shown in Fig. [TJ Here Er = is the recoil 
energy, where m is the boson mass. 

At zero temperature the many-body quantum state 
of interacting bosons in an optical lattice is well de- 
scribed by the Bose-Hubbard model when the lattice is 
so deep that only the lowest energy level of each lat- 
tice site is occupied and tunneling occurs only between 
nearest-neighboring sites [ll|, |l2j. To study bosons in ID 
DWOLs, we use the two-legged Bose-Hubbard Hamilto- 
nian (BHH) 

H = J2\ (^i^j.v - !) + ( n f ~ ^j,v 

j S^IL.i?.} 

- t \\( a ]+i,v a i,n+ h - c -)) - t±{a], R aj,L +h.c.) 



(2) 



,t 



a'j creates a boson at the lowest level localized on the left 



(right) of the j-th double-well for 77 = L(R), and fij jV = 
Oj „a,j,r) is the number operator. The lattice parameters 
can be related to the DWOL depths V[ n , V out , and V z and 
the recoil energy E R by assuming that 9 <§; 1 and that 
each well is sufficiently deep to approximate the Wannier 
function by a Gaussian. In this case, the on-site interac- 



tion is U ~ E R (8n) 
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offset) of the double wells is A ~ 
Here, a = V out /V in , s = V in /E R , 
is the s-wave scattering length, and 
f TZ^J ■ O is the curvature of the 
and /i is the chemical potential, which 

Notice that U, 
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parabolic trap, 

controls the total number of bosons N. 
tii, and t± are approximately independent of the tilt as 
long as 9 <C 1 ■ 

We now describe briefly the zero-temperature phase 
diagrams of the two-legged Bose-Hubbard model with no 
confining potential (Q = 0) @,[T3,[i3. In Figs. Ufa) -(c), 
we show the phase diagrams in the (fi/U, tn/U) plane for 
t±/U = 0.1. There we see the usual MI phase with unit 
filling, v = 1 (corresponding to two atoms per double 
well), and also one with half filling, v — 1/2 [a. IT3|. The 
SF-MI phase boundary in the (/j,/U,t\\/U) plane signifi- 
cantly depends on t±/U [H| and X/U 0. In particular, 
the unit-filling MI phase depends nonmonotonically upon 
X/U. As seen in the upper lobes of the phase diagrams, 
the v — 1 MI region shrinks initially as X/U is increased 
from zero. The size of the MI region is minimized at 
X = U, where the local states |1, 1) and |0, 2) are nearly 
degenerate; |j"Il,tir) designates the local Fock state with 
nL(nn) bosons on the left (right) of a double well. Due 
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FIG. 2: (color online) Results at three values of X/U: (a) 
and (d), 1 (b) and (e), and 2 (c) and (f). In the left three 
figures, we show the phase diagrams of the homogeneous sys- 
tem in the (fi/U,t»/U) plane for t± = 0.1(7, which have been 
calculated by means of the infinite-size version of the TEBD 
method. In the right three figures, the local number of bosons 
at the j-th double- well fij (squares), that at the j-th left well 
(diamonds), that at the j-th right well (triangles), and the 
fluctuation Oj of fij (stars) are shown. The dashed lines in 
the phase diagrams represent the values of the local chemical 
potential and the intrachain hopping realized in the right fig- 
ures. The upper (lower) end of the dashed lines corresponds 
to the center (edges) of the trapped gas. 



to this degeneracy, coherence is developed within each 
double well and the system favors the SF phase. As X/U 
is increased further, the MI region grows again. This 
nonmonotonic behavior leads to a reentrant phase tran- 
sition from a MI to a SF and again to a MI, induced by 
increasing X/U. 

Since we are mainly interested in the reentrant phase 
transition, we investigate how the ground state properties 
change with increasing X/U for the following fixed values 
of the other parameters; N = 98, fl/U = 0.001, t±/U = 
0.1, and t\\/U = 0.05. This value of t\\/U is sufficiently 
small such that the MI domain with v = 1 is present 
when A = or A 3> U, while it is so large that the v = 1 
MI domain is absent when A = U. The number of double 
wells is chosen as L = 85. This value of L is so large that 
the bosons do not see the edge of the system. 

To analyze the two-legged BHH (J2J , we use the finite- 
size version of the TEBD method [9], which provides 
a precise ground state for ID quantum lattice systems 
via propagation in imaginary time. While the maximum 
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FIG. 3: (color online) Density profile rij as a function of X/U. 
All the other parameters are fixed as N — 98, L = 85, tu/U = 
0.05, t±/U = 0.1, and Q/U = 0.001. 



number of bosons per site is N max = oo, convergence 
is already achieved in our numerical calculations when 
Amax = 5. We calculate the ground states of the Hamil- 
tonian ([2]) with x = 80, where x is the size of the basis 
set retained in the TEBD procedure Q. 

In Fig. 02 we show the density profile fij — rij^L+nj.R = 
^] (ftj^), namely the local number of bosons in the j-th 
double well as a function of X/U. For A = 0, the sys- 
tem is essentially a plateau with unit filling (fij = 2), 
with small SF regions with incommensurate filling at the 
edges. This plateau indicates the presence of an incom- 
pressible MI domain [1, M, EE EE HI- As X/U is in- 
creased, the unit-filling plateau starts to melt at a cer- 
tain point and a new MI plateau with half filling (fij = 1) 
emerges. In the vicinity of A = U, the unit-filling plateau 
has been completely melted away and the density profile 
for fij > 1 is smooth, with the reflected parabolic shape of 
the confining potential that is characteristic of SF phases 
in the regime where the local density approximation is 
valid [15(. As X/U is increased further, the unit-filling 
plateau forms again at the center of the system. Thus, 
as in the homogeneous case, the reentrant phase tran- 
sition is caused by increasing the tilt parameter in the 
presence of a parabolic potential [l9[ . 

For a better understanding of the reentrant phase tran- 
sition, three slices from Fig. are shown in Figs. [21(d)- 
(f) together with ri^z,, rij^R, and the fluctuation of fij, 

a j = yj (fij) — fij, which is small in the MI regions and 

relatively large in the SF regions. One can roughly in- 
terpret the spatial dependence of fij by a local density 
approximation. Introducing the effective local chemi- 
cal potential /x° ff = fx — Qj 2 , fij is approximated by 

n(/^ ff ), where n((J>) is the number of bosons per double 
well for uniform systems. Scanning the effective chem- 
ical potential along the dashed lines in the phase dia- 
grams of Figs. H{a)-(c), the behavior of the local density 
in Figs. [2jd)-(f) can be explained. When A — U, for ex- 
ample, the upper end of the dashed line in Figs. [U^b) is 
located in a SF phase with n > 2 as realized at the center 
of the trapped gas in Fig. [5£e) . As fi is decreased from 
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FIG. 4: (color online) (a) Quasi-momentum distribution S(k) 
as a function of X/U. (b) Visibility V (squares) and center 
density n J= o (triangles) as functions of X/U. 



the upper end, the dashed line enters the MI phase with 
half filling, corresponding to the half-filling plateau of the 
trapped gas. As /i is decreased further to the lower end, 
the dashed line exits to a SF phase with ri < 1, which is 
also realized at the boundary of the trapped gas. Notice 
that the local density approximation breaks down in the 
vicinity of the interfaces between the commensurate and 
incommensurate regions, where the singularities due to 
quantum critical behavior observed in uniform systems 
are removed [ljj • 

In Figs. [2d)-(f), the population of bosons in the left 
wells shifts gradually to the right wells as X/U is in- 
creased. At A = U, almost all the bosons occupy the 
right wells, namely rij^L — 0, for fij < 1. In contrast, a 
considerable amount of bosons remains in the left well as 
rij.L — " J 2 1 for fij > 1, since the local state with unit- 
filling is approximately given by (|1, 1) + |0, 2))/^/2 due to 
the competition between the onsite interaction and the 
tilt. At A = 2U , the population of bosons in the left well 
almost vanishes, namely rij^ — 0, in the entire system. 

Next, we calculate the quasi-momentum distribution 
S(k) = L- 1 J2ji E„ e ifc tf-^(at n oi )1? ), where d is the lat- 
tice spacing. Since the true-momentum distribution is 
expressed as the product of |w(fc)| 2 and S(k) (l8j |. where 
w(k) is the Fourier transform of the Wannier function 
in the lowest Bloch band, the quasi-momentum distri- 
bution can be extracted by dividing the true-momentum 
distribution observed in experiments by |w(fc)| 2 @. S(k) 
is shown as a function of X/U in Fig. SJa), where the 
k = peak is the sharpest in the vicinity of A = U. 
This means the intrachain coherence increases as the SF 
region around the trap center emerges. 

To quantify the sharpness of the peak in S(k), we cal- 
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FIG. 5: (color online) (a) Local coherence between the left and 
right wells (a] L aj,R) for X/U = (squares), 1 (diamonds), 
and 2 (triangles), (b) Total coherence C versus X/U . 

culate the visibility of interference patterns defined as 
V = (S max - 5 mi „)/(5 max + S min ) PHI. In Pig. Hp)), 
we show V and the center density nj =0 as a function of 
X/U. As indicated by the first arrow, the visibility ex- 
hibits a sudden jump at the point where the unit-filling 
MI plateau starts to melt. In the vicinity of A = U, where 
the unit-filling MI plateau does not exist, the visibility 
is distinctively large compared to the regions where the 
unit-filling MI plateau is present. As X/U is increased 
further, the visibility decreases the most rapidly at the 
point indicated by the second arrow, where the v = 1 
MI plateau forms again. Notice that the change of the 
visibility at the second arrow is not as distinctive as the 
jump at the first arrow because the v = 1 Mott plateau 
at A = 2Z7 is small compared to that at A = 0. This ten- 
dency of the visibility suggests that the reentrant phase 
transition can be characterized in experiments through 
the observation of matter-wave interference patterns. 
Another key observable in DWOLs is the coherence 

between the left and right wells C = X^ ( a j L a j,RI- The 
coherence C can be extracted from the double-slit diffrac- 
tion pattern, which has been observed in the time-of- 
flight images in the horizontal xy plane [j], 0, 3- The 
larger C is, the clearer the diffraction pattern is. 



We calculate first the local coherence (aj L aj.B.) as 
shown in Fig. Efa). At A = 0, since the local state with 
half-filling is well approximated by the bonding state, the 
SF regions at the boundary of the trapped gas have large 
local coherence. As X/U is increased, the region with 
fij < 1 immediately loses the local coherence. In con- 
trast, in the region with fij > 1 the local coherence is 
pronouncedly large at A = U . The development of lo- 
cal coherence is due to the degeneracy of |1, 1) and |0, 2) 
states, and it drives the system into the SF phase. As 
X/U is increased further, the entire system tends to lose 
coherence. In Fig. E^b), the total coherence C is shown 
as a function of X/U . While the total coherence is sig- 
nificantly reduced when A is increased from U , it hardly 
changes in the region A < U. This happens because the 
gain of local coherence for fij > 1 and the loss for fij < 1 
almost cancel each other out. Consequently, the double- 
slit diffraction pattern can not capture the signature of 
the reentrant phase transition. 

In conclusion, we have studied the SF-MI transition 
of parabolically trapped Bose gases in a ID double-well 
optical lattice by using the time-evolving block decima- 
tion (TEBD) method. We have calculated the density 
profile as a function of the tilt of the double wells and 
shown that a reentrant phase transition between MI and 
SF phases occurs in the presence of a parabolic confine- 
ment. We have calculated also the quasi-momentum dis- 
tribution and found that the matter-wave interference 
pattern, which is one of the most common observables in 
experiments with cold atoms, contains sufficient informa- 
tion to characterize the reentrant phase transition. Wc 
would like to emphasize that, unlike the results of mean- 
field theories, our results based on the TEBD method are 
quantitative, thus we expect the reentrant phase transi- 
tion to be observed in future experiments near the pa- 
rameter values discussed in this paper. 
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